Localization of Two Interacting Particles in One-Dimensional Random Potential 
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We investigate the localization of two interacting particles in one-dimensional random potential. 
Our definition of the two-particle localization length, ^, is the same as that of v. Oppen et al. [Phys. 
Rev. Lett. 76, 491 (1996)] and ^'s for chains of finite lengths are calculated numerically using the 
recursive Green's function method for several values of the strength of the disorder, W, and the 
strength of interaction, U. When (7 = 0,^ approaches a value larger than half the single-particle 
localization length as the system size tends to infinity and behaves as ^ ~ W~'^° for small W with 
— 2.1 ±0.1. When U 7^ 0, we use the finite size scaling ansatz and find the relation ^ ~ W"" with 
= 2.9 ± 0.2. Moreover, data show the scaling behavior ^ ~ W'"" giblUl/W^) with A = 4.0 ± 0.5. 
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PACS number(s); 72.15.Rn, 71.30.+h 

Recently, there has been intensive attention [0-0 fo- 
cused on the problem of the localization of two inter- 
acting particles in one-dimensional (ID) random poten- 
tial. With a few assumptions on the statistical nature 
of single-particle localized states, Shepelyansky |^ has 
mapped the problem approximately to a random band 
matrix model and obtained an expression for the two- 
particle localization length, ^, as 
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where U is the on-site interaction in unit of the hop- 
ping energy between nearest neighbor pair sites, and ^1 
the single-particle localization length. This expression is 
surprising because it implies that ^ can exceed ^1 at suf- 
ficiently small disorder, i.e. sufficiently large ^1. Later 
Imry B] has provided a support for Eq. (1) by invoking 
the Thouless scaling argument. However, the methods 
employed in Refs. [1] and [2] are partly approximate and 
the strict validity of the expression of Eq. (1) is question- 
able as discussed in^.g. Refs. [3-8] and [10]. 

Many authors [^-§1 have tried to find more refined ex- 
pressions than Eq. (1) by improving the assumptions of 
Shepelyansky. However, at this stage, there exist con- 
troversies yet as to the quantitative expression for ^ like 
Eq. (1). Frahm et al. obtained the relation ^ ^^-^^ 
by the transfer matrix method while an approximate cal- 
culation of Green function by v. Oppen et al. [^ leads 
to the hypothesis = ^i/2 + c|[/|^2: where c is a con- 
stant depending on the statistics of the particles. With 
the assumption that the level statistics of two interact- 
ing particles is described by a Gaussian matrix ensemble, 
Weinmann and Pichard Q argued that ^ increases ini- 
tially as |t7| before eventually behaving as U"^. Moreover, 
very recently, Romer and Schreiber have claimed the dis- 
appearance of the enhancement as the system size grows 
(see Refs. [7] and [8]). 



Some of these discrepancies, especially between nu- 
merical studies, are due to different definitions for two- 
particle localization length between authors and also to 
lack of careful analysis of the finite size effect of the sys- 
tem size. The system under study is a "quantum me- 
chanical two-body problem" in a sense. Motion of the 
two particles can be decomposed into the motion of the 
center of mass (CM) and that of the relative coordinate. 
We are interested in the CM motion since the wavefunc- 
tion describing the relative motion would not be differ- 
ent from that arising from the single-particle localization 
problem in the thermodynamic limit if the interaction is 
short-ranged. Therefore, in this paper, we use the same 
definition for ^ as introduced by v. Oppen et al. for 
the measure for localization length of the CM: 
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Here, G is the Green function and \i,j) is a two-particle 
state in which the particle 1 (2) is localized at a site 
i (j). The above definition is reasonable for a descrip- 
tion of the CM motion as long as U is smaller than or 
of the order of the hopping energy between sites [1). 
In practice, we calculate £,n defined below in Eq. (4) 
for chains of finite lengths without any approximation 
for several values of W and U. We then estimate ^ 
by extrapolating using the finite size scaling ansatz. 
When U = 0, we find ^ - W-"" with i^q = 2.1 ± 0.1. 
Data for U ^ lead to the relation ^ ^ W~'^ with 
ly = 2.9 ± 0.2. Also the data lead us to propose a scal- 
ing form ^ ~ W~''°g{b\U\/W^), where g{y) is a scaling 
function with the property g{y — s- 0) = constant and 
g{y ^ 00) ~ yi^-i^o)/^^ A is given as 4.0 ± 0.5. 

We work within the tight-binding equation given by 
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where ipi. 



E is the energy of the two par- 



ticles and (5m, n the Kronecker delta, m and n are the 
site indices of a chain of length N and range from 1 to 
N , em is the random site energy chosen from a box dis- 
tribution with interval [-W/2,W/2] ^ and the hard 
wall boundary condition, i.e. "00 n = and so on, is 
used. As was previously noted if one interprets 

(to, n) as Cartesian coordinates of a square lattice of 
size N X N, the Hamiltonian describes a single parti- 
cle in a two-dimensional random potential. In Eq. (2), 
the thermodynamic limit is first taken and then the limit 
|n — m| — > oo. To estimate this quantity, we define a 
sequence ^at as 

^ < ^T^ln|(l,l|Gjv|Af,iV>| >, (4) 
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where Gn represents the Green function for a chain of 
length N and the double brackets represent the con- 
figurational average. To be specific, calculation of Gjv 
amounts to evaluation of the inverse of the matrix, 
(E — Ti), the size of which is iV^ x N"^. One can calcu- 
late several elements of Gat, i.e. the elements involving 
the sites of two opposite edges of the square lattice, very 
efficiently using the recursive algorithm of MacKinnon 
and Kramer |12|. We assume that approaches ^ as 
N ^ oo. 

The on-site interaction of the Hamiltonian given by 
Eq. (3) is relevant only to the spatially symmetric states, 
which would be realized, say, for a pair of electrons with 
total spin zero. One can easily see that the contributions 
to Eq. (2) are only from the spatially symmetric states 
from following consideration. The Green's function rep- 
resents the transition amplitude from an initial state to a 
final state and since the Hamiltonian, Eq. (3), is invari- 
ant under the exchange operation of two particles, the 
parity of the wavefunction is conserved during the time 
evolution. Since the initial state of Eq. (2) is a doubly 
occupied state, i.e. a spatially symmetric state, we are 
treating only the contributions from symmetric states. 

Numerical calculations of £,n for various values of W, U 
and N are performed for E — Q without approximation. 
N is varied within the range 10 < < 200 and for a 
given parameter set, configurational average is performed 
over sufficiently many different realizations to control the 
uncertainties of within 1%. 

We first examine the case oiU = 0, i.e. the noninter- 
acting two particles. In this case, when the total energy 
of the system is fixed to E^ the two-particle wavefunction 
is a superposition of the products of two single-particle 
states of energy E' and E — E', and the Green function 
is given by the convolution of two single-particle Green 
functions as 

{z,i\G{E)\j,j) ^ I dE'{i\G^{E')\j){i\G^{E-E')\j). 



(5) 

It is a nontrivial problem to calculate <^(J7 — 0) since 
there exist contributions from various energies. Some 



authors |Jj5|] have assumed the relation ^([/ = 0) = Ci/2, 
i.e. half the single-particle localization length, which 
should be, however, seriously examined. Our numerical 
data presented in Fig. 1 show that the assumption is not 
strictly valid. The filled symbols on the N — oo axis rep- 
resent Ci/2 calculated from the expression ~ 105/1^^ 
p3[ while the empty symbols are our numerical results 
for ■ Taking into account the fact that the uncertainty 
of each data point is smaller than the symbol size, ^a? 
does not seem to extrapolate to ^i/2 as N tends to in- 
finity. Moreover, the discrepancy between the two quan- 
tities becomes larger as W gets smaller. Therefore, we 
conclude that within the definition of Eq. (2), the single- 
particle localization length is not an adequate parameter, 
if it is qualitatively, for a quantitative description of two- 
particle localization problem. From the data of = 200, 
we get ^([7 = 0) ~ 70/11^''° with = 2.1 ± 0.1. 

Next, we discuss the case oiU 7^ 0. Figure 2(a) shows 
the results for U = 1.0 and W ranging from 0.5 to 10.0. 
The y axis label represents the renormalized localization 
length, i.e. ^at divided by the system size. For larger val- 
ues of W and A, £,n/N behaves as ^ l/N, which implies 
the convergence of ^Af's to their constant limiting values. 
This means that the condition A^ ^ ^ is well satisfied 
for these data. However, for smaller values of W, i.e. for 
W ranging from 0.5 to 1.5, it is not easy to deduce the 
value of ^ since S^n 's increase steadily within the range of 
A^ presented. Therefore we rely on the scaling idea [ p2[ , 
which states that £,n /A is given by a function of a single 
parameter, i.e. N/^: 



^N/N = fiN/0. 



(6) 



The implication of Eq. (6) is that on a log-log plot all 
data points of Fig. 2(a) fall on a single curve when trans- 
lated by In^(M^) along the x axis. As a result, ^{W)'s 
can be obtained as fitting parameters. The result of data 
collapsing is shown in Fig. 2(b) for the data set A > 50. 
^{W = 5.0) has been obtained to be 2.87 ± 0.01 by fit- 
ting the data set for W — 5.0 and A^ > 50 to the formula 
= £. — A/N iQ, where A is a constant. Other re- 
maining values of ^(W)'s are obtained by examining the 
amount of relative translations with respect to the data 
set of W = 5.0. The scaling plot is quite good and one 
can see that the scaling function f{x) behaves as 



1/^ 
1/x 



if X < 1, 
if a; > 1. 
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As was previously mentioned, the asymptotic behavior 
for a; 3> 1 represents the convergence of £,n's to On 
the other hand, the behavior for a; <C 1 is very interesting 
since the same asymptotic behavior has been found for 
noninteracting disordered ID systems |13|. For the non- 
interacting case, the resistance of a chain of length A 
is related to the single-particle localization length as 
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= [cosh(2Ar/^o,) _ 1^/2. (8) 

For N/S^^ ^ 1, the right hand side of Eq. (8) reduces 
to ~ (N/S.n)'^ ^ ^/Ci- Therefore for the noninteract- 
ing case, the asymptotic behavior for TV <C represents 
the metalHc behavior of the resistance, i.e. the hnear in- 
crease of the resistance as the chain length in the metal- 
he regime. Though no explicit expression like Eq. (8) is 
not present for the system under study in this paper, we 
believe that the same asymptotic behavior for the two 
cases found here is a strong indication that the definition 
in Eq. (2) is a physically reasonable one. 

The ^'s thus obtained as a function of W are plotted 
in Fig. 2(c) For 0.5 < < 5.0 they are reasonably 
well fitted to the form of ~ W~'^, where ly is given by 
2.9 ± 0.1. Within the error, this value for i/ is different 
from 1^0 = 2.1 ± 0.1, i.e. the critical exponent for U — 0, 
and from 4.0, which is the value expected by Eq. (1). 

Further calculations and similar scaling analyses have 
been performed for other values of U, i.e. 0.2, 0.5, 0.7 
and 1.5 upto system size N — 200. It is difficult to deter- 
mine ^ for W < 1.5 and U < 1.0 since the corresponding 
data of ^Tv's do not show scaling behaviors due to the 
fact that sufficiently large system sizes have not been 
reached for these parameters. The resulting ^'s (for 1.5 
< < 5.0 if [/ < 1.0 and for 0.7 <W < 5.0 if f7 > 1.0) 
give V = 2.7, 3.0, 2.9 and 3.1 for U = 0.2, 0.5, 0.7 and 1.5, 
respectively. Since we do not expect that v depends on 
U, we interpret the variations of the values for v as re- 
sulting from numerical uncertainties. Therefore our final 
result for the critical exponent of ^ is 2.9 ± 0.2. 

Our result for v implies that introduction of nonzero 
U changes the critical behavior of ^ and, in analogy with 
thermal critical phenomena, the point W — U — may 
be regarded as a multicritical point and the line W = 
as a critical line in the W — U plane. Then, one may 
assume a scaling form for ^ as follows; 

^^W-''"gib\U\/W^), (9) 

where g{y) is a scaling function, A a crossover exponent 
and b a constant. Here, we used the fact that Eq. (3) 
is symmetric for £^ = so that ^ depends only on the 
absolute value of U. The scaling function should satisfy 
g{y ^ 0) = constant and g{y oo) ~ y('^-^o)/A j^j. 
sistency. We obtain reasonably good scaling plots within 
the range A = 4.0 ± 0.5. The scaling plot for A = 4.0 
is shown in Fig. 3, where (,W^'^ is plotted against U /W^ 
for various values of W and U . Although the data for 
W < 1.5 may appear to deviate from the scaling curve, 
taking into account rather large numerical uncertainties 
of these data, one can expect that they are consistent 
with the scaling behavior of other data points. We ex- 
pect that the crossover between the two asymptotic be- 
haviors occurs at ?/ ~ 1 so that the constant b is estimated 
to be of the order of 100 from Fig. 3. Data within the 
range 0.01 < U/W^ < 1.0 (1 < y < 100) approximately 
obey the form ~ y'^'^^, which is shown as a straight line. 



Since we expect the asymptotic behavior ^ yi'^-'^o)/A -f^j, 
this regime, we obtain {ly — i'o)/A ~ 0.23, i.e. ly ~ 3.0, 
which is in good agreement with our previous estimate, 
i.e. ly = 2.9±0.2. At this stage, we have not found a phys- 
ical mechanism regarding the scaling parameter, U/W^ 
with A = 4.0 ± 0.5. The quantity ^{U) - £,{U = 0) might 
be also of interest. Our data show that this is consistent 
with a form ^{U) - ({U = 0) - W-^-^iU/W^y^^ in the 
region bU/W^ < 1.0. This behavior shows that the first 
correction to g{y) for y ^ 1 is given as ^ y^^^. On the 
other hand, our result confirms the enhancement of two 
localization length due to the interaction. The arrow of 
Fig. 3 represents the value of ^W^"" for ?7 = so that it is 
clearly seen that ^ monotonically increases with respect 
to U. 

Finally, we point out differences between our work and 
some of those previously reported. References [3] and [7] 
deal with exactly the same system as ours but use a differ- 
ent definition for the two-particle localization length. As 
mentioned before, the problem can be considered as that 
of a noninteracting single particle in a two-dimensional 
potential. These authors study the evolution of a state 
along one of the edges of the square lattice. However, 
in this paper, we are concerned with the "pair propaga- 
tion" , more generally, the propagation of the CM of the 
two particles. The definitions for ^ given by Eqs. (2) and 
(4) describes the propagation along one of the diagonals 
of the square lattice, instead of that along the edge. Our 
definition of ^ is exactly the same as that of v. Oppen 
et al. 1^. However it should be noted that in their work, 
calculation of ^ involves an approximation; the approxi- 
mation scheme used in Ref. [5] fails for small values of U 
while our results are valid for all values of U . 

In summary, we have investigated numerically the lo- 
calization of two interacting particles in ID random po- 
tential using the definition introduced previously for the 
two-particle localization length. While we find the en- 
hancement of ^ by the interaction, critical properties of ^ 
are different from those reported in previous studies. We 
ascribe the differences to the approximation used in one 
case and, in the other cases, to different definition of ^. 
Further works are needed to connect the resistance and 
the two-particle localization length and to elucidate the 
relation between ^ and . 
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FIG. 1. ^iv (open symbols) as a function of 1/A^ for 
U — and ^i/2 (filled symbols) calculated from the expression 
^1 ~ W5/W'^: W = 1.0 (box), 1.5 (circle), 2.0 (uptriaiigle), 
3.0 (downtriangle), 4.0 (diamond) and 5.0 (pentagon), from 
top to bottom. The uncertainty of each data point is less than 
the symbol size. 



1000 



100 

U = 



U = 0.2 o 

0.5 ■ 

0.7 o 

1.0 • 

1.5 ^ 



10 

0.0001 0.001 0.01 




0.1 



10 



100 



FIG. 2. (a) £,n/N as a function of N for U = 1.0. (b) 
Scaling plot constructed from the data of (a) for N > 50. 
(,{W)'s are obtained as fitting parameters by this procedure, 
(c) ^ as a function of W. The data for W < 5.0 fit well to a 
straight line of ~ W~^'^ as shown. 



FIG. 3. The scaling plot of ^W" versus U/W^ with 
vo = 2.1 and A — 4.0. The typical uncertainty of data for 
U/W'^ > 1.0 {W < 1.5) is shown at the rightmost data point, 
while that for U/W^ < 1.0 {W > 1.5) is shown at the left- 
most data point. The straight line is ~ x"'^^ and the arrow 
represents ~ 70 for U = 0. 
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